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Abstract. Electron heating and ionization dynamics in capacitively coupled radio 
frequency (RF) atmospheric pressure microplasmas operated in helium are investigated 
by Particle in Cell simulations and semi-analytical modeling. A strong heating of 
electrons and ionization in the plasma bulk due to high bulk electric fields are observed 
at distinct times within the RF period. Based on the model the electric field is 
identified to be a drift field caused by a low electrical conductivity due to the high 
electron-neutral collision frequency at atmospheric pressure. Thus, the ionization 
is mainly caused by ohmic heating in this "Sl-mode". The phase of strongest bulk 
electric field and ionization is affected by the driving voltage amplitude. At high 
amplitudes, the plasma density is high, so that the sheath impedance is comparable 
to the bulk resistance. Thus, voltage and current are about 45° out of phase and 
maximum ionization is observed during sheath expansion with local maxima at the 
sheath edges. At low driving voltages, the plasma density is low and the discharge 
becomes more resistive resulting in a smaller phase shift of about 4°. Thus, maximum 
ionization occurs later within the RF period with a maximum in the discharge center. 
Significant analogies to electronegative low pressure macroscopic discharges operated 
in the Drift- Ambipolar mode are found, where similar mechanisms induced by a high 
electronegativity instead of a high collision frequency have been identified. 



PACS numbers: 52.20.-j, 52.25.Jm, 52.27.Cm, 52.50.-b, 52.65.Rr, 52.80.Pi 
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1. Introduction 

Microscopic capacitively coupled plasmas operated at atmospheric pressure and at a 
radio frequency of typically 13.56 MHz, e.g. microscopic atmospheric pressure plasma 
jets (/x-APPJ), are frequently used for surface processing and medical applications [TH5]. 
The non-thermal glow discharge plasma in these sources consists of hot electrons and 
cold heavy particles (ions, neutrals) close to room temperature [6]. Such discharges 
can be used for sensitive surface treatments including human tissue. Atmospheric 
pressure plasmas avoid the necessity of expensive vacuum systems required for processing 
applications in low pressure discharges, while providing a high degree of dissociation and 
an effective generation of reactive species useful for surface treatment [7HT2]. 

There are different types of APP Js such as coaxial and plane parallel configurations [2||5] . 
The latter design concept is based on the plasma jet introduced by Selwyn et al. in 
1998 [13] and modified by Schulz-von der Gathen et al. [H]: The feed gas flows between 
two electrodes separated by a gap of about 1-2 mm and driven at 13.56 MHz. In the 
experiment the electrodes are typically made of stainless steel and enclosed by quartz 
windows including the plasma volume and the effluent. In this way direct optical access 
to the plasma and the effluent is provided. Usually the discharge is operated in helium 
with some optional admixture of oxygen and/or nitrogen with typical gas velocities 
around 100 m/s [T5] . 

The generation of reactive species, that determine surface processes in the effluent, 
is caused by electron impact excitation, ionization, and dissociation in the plasma 
volume between the electrodes. Thus, a detailed understanding of the dynamics of 
highly energetic electrons in the plasma is essential and provides the basis for any 
optimization of surface processing applications. The electron dynamics in /z-APPJs has 
been investigated experimentally by, e.g. Schulz-von der Gathen et al. [T4P15]. Benedikt 
et al. [16], and Kong et al. [17] by Phase Resolved Optical Emission Spectroscopy [18] 
as well as numerically by Waskoenig et al. [19] and Kong et al. [20Tl2~2"] . Within the 
RF period several maxima of the emission at distinct positions and times have been 
observed: (i) During sheath expansion electrons are accelerated towards the opposing 
electrode and cause excitation/ionization adjacent to the expanding sheath edge [22J. 
(ii) At the time of maximum sheath voltage excitation and ionization by secondary 
electrons is observed at both sheaths [I9l[23l|24j. (iii) During sheath collapse another 
excitation maximum is observed in /i-APPJs at each electrode [T7J. This maximum 
is assumed to be caused by an electric field reversal localized at the sheath edge and 
caused by electron-neutral collisions [25]. (iv) Significant excitation and ionization inside 
the bulk at the times of fastest sheath expansion are observed [T6"|ITT]. These maxima 
have been correlated with high bulk electric fields, but their exact physical origin is 
not completely understood. At low driving voltage or power, the excitation during 
sheath expansion is typically stronger than the excitation by secondary electrons [17J. 
Such mode transitions induced by changing the RF voltage amplitude are similar to the 
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mode transitions discussed by Belenguer and Boeuf at low pressures [25] . 

Here, we investigate the electron heating and ionization dynamics in an atmospheric 
pressure microplasma with plane parallel electrodes driven at 13.56 MHz as a function 
of the RF voltage amplitude by Particle in Cell (PIC) simulations and semi-analytical 
modeling. The discharge is operated in helium. We reveal the origin of the ionization in 
the bulk and show, that the ionization maxima adjacent to the collapsing sheaths are not 
caused by a classical localized field reversal under the conditions investigated. Based on 
the analytical model, we demonstrate, that the strong ionization in the bulk and at the 
sheath edges is caused by a high electric field inside the bulk at the time of maximum 
current. This high field originates from a low DC conductivity due to a high electron- 
neutral collision frequency at atmospheric pressure in the bulk. The phase shift between 
current and voltage is found to be affected by the driving voltage amplitude. Thus, 
maximum electron heating and ionization occur at different times within the RF period 
depending on the voltage amplitude, that also affects the spatial profile of the electron 
heating and ionization rates. We compare our results to low pressure macroscopic 
electronegative capacitive discharges operated in CF 4 , where similar effects are caused 
by the high electronegativity instead of a high collision frequency [27]. We conclude, 
that a novel mode of discharge operation, the f2-mode, is present in atmospheric pressure 
microplasmas, where ionization is dominated by ohmic heating in the bulk. 

The paper is structured as follows: In section 2, the PIC simulation and the semi- 
analytical model to describe the electric field in the bulk are introduced. In the third 
section, the results are presented and discussed. Finally, conclusions are drawn in section 
4. 

2. Methods 

2.1. Particle- In- Cell Simulations 

We use a ld3v (one spatial dimension and three dimensions in velocity space) explicit 
kinetic code based on the particle-in-cell algorithm to study the heating and ionization 
dynamics in capacitive atmospheric pressure microplasmas with plane parallel electrodes 
of identical surface areas separated by a gap of 2 mm. Collisions are treated by the 
Monte-Carlo method. The code works in the electrostatic approximation, i.e. E = — V0. 
The driving voltage waveform is = <PqCos{2'k ft) with f = 13.56 MHz and 0o — 500 
V, 330 V. To overcome the limitation of a very small time step - the electron elastic 
collision frequency has to be resolved - the simulations are accelerated: We use a coarse- 
sorting algorithm for massive parallelization of the code on graphics processing units 
(GPU) [28J. We consider electrons and He + ions taking cross sections from [291130] 
and use an ion- induced secondary electron emission coefficient of 7 = 0.1 as well as an 
electron sticking coefficient at the electrodes of 0.5. The gas temperature is set to 350 
K. 
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In order to compare the results obtained for atmospheric pressure microplasmas operated 
in helium to low pressure macroscopic electropositive and electronegative discharges 
we also perform PIC simulations of Ar and CF4 discharges with an electrode gap of 
1.5 cm operated at 80 Pa. The cross sections for Ar are taken from [30]. In case 
of CF 4 , we consider electrons and the ions CF^ , CFJ, F~ using cross sections and 
rate coefficients from [3T1[32]. A rate constant of 5.5xl0 _13 m 3 s _1 is used for the ion- 
ion recombination (CF^ + CF3 , CF3" + F~). The coefficient for secondary electron 
emission at the electrodes due to ion bombardment is varied in the electropositive case 
and set to 7 = 0.1 in electronegative CF 4 . The probability of sticking of electrons at the 
electrodes is 0.8. In both the electropositive and electronegative cases at low pressure 
the gas temperature is kept constant at 350 K. [33]. 



2.2. Semi- analytical model for the bulk electric field 

In order to understand the physical origin of the bulk electric field obtained from the 
PIC simulations, an analytical expression for the electric field is deduced, following [25J . 
The model is based on a combination of the electron continuity and momentum balance 
equations, i.e.: 

**" dT ' + S„ (1) 



dt dx 
8T e d(v e T e ) 1 dp e e 

-sr = ~ -k n e E - u ce T e . (2) 

dt dx m e dx m e 

Here, p e is the electron partial pressure, with p e = n e T e , n e is the electron density, T e 
is the electron temperature, E is the electric field, r e is the electron flux, m e is the 
electron mass, S e is the ionization source, and v ce is the frequency of elastic electron- 
neutral collisions. Combining ([Q) and ([2]) and assuming quasineutrality in the bulk as 
well as a spatially homogeneous electron temperature yields: 

m c (dj e . 1 / , 2 , 2 x dn e j e S e \ 



Here, j e = en e v e is the electron conduction current density and j t h, e = en e yT e /m e . The 
first and third term of equation [3] correspond to electric fields caused by inertia effects, 
while the second term is a drift field, the fourth term is the ambipolar field, and the 
fifth term corresponds to an electric field caused by an ionization source. 



3. Results 



Figure [T] shows different plasma parameters resulting from PIC simulations of 
capacitively coupled microscopic atmospheric pressure plasmas operated in helium with 
an electrode gap of 2 mm at 13.56 MHz and a voltage amplitude of 0o — 500 V. The 
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Figure 1. PIC simulation results (He, 1 atm, 13.56 MHz, 2 mm gap, 500 V): Spatio- 
temporal plots of the (a) electron density, (b) space charge density, (c) electric field, 
(d) electron conduction current density, (e) electron heating, and (f) ionization rate, 
(g) shows the electric field profile at the time of max. ionization [vertical dashed lines 
in (a) - (f)] and the second term (drift field) of eq. ((3]). (h) shows the conduction 
current density in the discharge center and the applied voltage as a function of time. 
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spatio-temporal distributions of the electron density, space charge density, electric field, 
electron conduction current density, electron heating rate, and ionization rate obtained 
from the simulation are shown in plots (a) - (f), respectively, while the electric field 
profile at the time of maximum ionization marked by vertical dashed lines in (a) - (f) 
resulting from the simulation and the analytical model [equation fl3])] is shown in plot 
(g). The electron conduction current density in the discharge center and the applied 
driving voltage waveform are shown in plot (h). 

electron heating rate ionization rate electric field electron density 
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Figure 2. PIC simulation results: spatio-temporal plots of the electron heating 
rate (first column), ionization rate (second column), electric field (third column), and 
electron density (fourth column) in Ar and CF4 discharges driven at 13.56 MHz and 
80 Pa with an electrode gap of 1.5 cm. First row: Ar, 100 V, 7 = 0. Second row: Ar, 
200 V, 7 = 0.2. Third row: CF 4 , 400 V, 7 = 0.1. The color scales are given in units of 
10 5 W m~ 3 (heating rate), 10 21 m~ 3 s" 1 (ionization rate), 10 3 V m _1 (electric field), 
and 10 15 m~ 3 (electron density). Reproduced from Ref. [27] . 

Analogous simulation results for the spatio-temporal electron heating rate, ionization 
rate, electric field, and electron density obtained in macroscopic low pressure capacitively 
coupled RF (CCRF) discharges operated in argon and CF4 at 13.56 MHz and 80 Pa 
with an electrode gap of 1.5 cm are shown in figure |2j The first row shows results for a 
discharge operated in argon at 80 Pa, 100 V voltage amplitude, and 7 = 0, the second 
row corresponds to an argon plasma at 80 Pa, 200 V, and 7 = 0.2, and the third row 
shows results for CF 4 at 80 Pa, 400 V, 7 = 0.1. The rows of figure |2] correspond to 
the 3 different modes of electron heating known in macroscopic low pressure capacitive 
RF plasmas: The first row corresponds to the a-mode, where electrons are heated by 
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Figure 3. Profiles of (a) the electric field obtained from the simulation (black solid 
line) and §S§ (dashed red line); (b) the second (drift field) and fourth (ambipolar held) 
terms of ([3]) in the bulk at t « 26 ns [vertical dashed lines in Fig. EJi)-(l)]- Discharge 
conditions: CF 4 , 13.56 MHz, 400 V, 80 Pa, 7 = 0.1. Reproduced from Ref. [27]. 



direct interaction with the expanding sheaths and cooled by direct interaction with 
the collapsing sheaths [3H436] . If the secondary yield, 7, is increased, the discharge 
is operated in 7-mode at otherwise similar conditions and the ionization is dominated 
by secondary electrons generated by ion impact at the electrodes. These 7-electrons 
are accelerated and multiplied effectively inside the sheaths at the times of maximum 
sheath voltage (second row, [38l[39]). Electronegative discharges can be operated in the 
Drift- Ambipolar (DA) mode, where ionization caused by strong electric fields inside the 
bulk at distinct times within the RF period dominates. This mode is observed in CF4 
discharges at 80 Pa such as shown in the third row |27j . 

In microscopic atmospheric pressure plasmas (figure [1]) , the electron density peaks in 
the central bulk region and decreases monotonically towards the electrodes similar to 
electropositive macroscopic CCRF discharges (1st and 2nd rows in figure [2]). In contrast 
to low pressure electropositive discharges the electric field is high inside the plasma bulk 
at two distinct times within the RF period. One phase of high bulk field is marked 
by a vertical dashed line in figure HJ The second phase of high bulk electric field 
occurs half an RF period later. At these times, the conduction current is high in 
the discharge center and maximum electron heating as well as ionization adjacent to 
the sheath edges are observed. Weak additional ionization by secondary electrons is 
observed inside the sheaths at times of maximum sheath voltage. This is similar to 
low pressure electronegative discharges, where also a high bulk electric field is found at 
distinct phases in the RF period (3rd row in figure EJ). 

Although the spatio-temporal ionization dynamics during sheath expansion in 
atmospheric pressure microdischarges [figure [1] (f)] looks similar to results obtained 
in low pressure macroscopic CCRF discharges operated in a- mode [figures [2] (b)], the 
physical mechanisms causing these maxima are different. At low pressures, electrons 
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directly interact with the expanding sheaths and are heated stochastically. Electron 
beams [35] are generated by the expanding sheaths and propagate into the bulk, where 
they cause ionization typically within a distance of one electron mean free path away 
from the sheath edge. This direct interaction of electrons with the time dependent 
sheath electric field results in cooling during sheath collapse. At atmospheric pressure, 
the electron mean free path is below 1 fim. However, ionization at the time of sheath 
expansion is observed up to 500 fim away from the expanding sheath edge. This 
ionization is caused by a high bulk electric field such as shown in fig. 1 (c), which 
has local extrema adjacent to the sheath edges at both electrodes. The semi- analytical 
model [equation (jHJ)] reproduces the electric field well at the time of maximum ionization 
using input parameters from the simulation in terms of j e , n e , T e , and v ce . By separating 
the contributions of the individual terms in equation ([3]), the physical origin of this high 
bulk field is revealed. It turns out, that the second term (drift field) dominates compared 
to the others. Thus, the comparison between model and simulation shows [figure [1] (g)], 
that the high bulk field is predominantly a drift field caused by a low conductivity 
due to the high electron- neutral collision frequency at atmospheric pressure. This field 
increases slightly from the discharge center towards the sheaths due to the decreasing 
plasma density at constant conduction current density. Any slope of the electric field 
profile is correlated with a small local space charge shown in figure [1] (b) jlQ] . Electrons 
are accelerated in this high bulk electric field and cause ionization, where the field 
is maximum, i.e. at the sheath edges. This also explains the local maximum of the 
ionization at the collapsing sheath edge, where no localized field reversal is present. 
Thus, the ohmic bulk heating of electrons in this fi-mode dominates the electron 
heating and ionization dynamics under these conditions. A remarkable difference to 
the electron heating dynamics in low pressure macroscopic CCRF discharges is the 
absence of electron cooling by the collapsing sheaths. This verifies the absence of direct 
interaction of electrons with the moving sheaths and the importance of the bulk heating. 

Under the conditions investigated here, current and voltage are 46.8° out of phase due 
to the high plasma density, which causes the sheath impedance to be comparable to the 
bulk resistance, i.e. 1/(u>rfC s ) ~ Rb- Here, ojrf = 27r • 13.56 MHz, C s ~ e^A/s is the 
effective capacitance of both sheaths determined from the time averaged sheath width, 
s, and Rt = (u ce db) / (cjp e e A) is the bulk resistance determined from the the bulk length, 
db, and the electron plasma frequency, u pe [37J. A is the electrode area. 

The presence of a high bulk electric field in electropositive atmospheric pressure 
microdischarges is similar to the situation in low pressure macroscopic and 
electronegative CCRF discharges operated in the Drift-Ambipolar mode (3rd row in 
figure E]), where a low conductivity is caused by a low electron density due to the high 
electronegativity [27]. Under these conditions the second and fourth terms of equation 
((3]) reproduce the high bulk electric field such as shown in figure |3j The shape of the 
ambipolar field is caused by the presence of local maxima of the electron density in the 
electropositive edge regions of this electronegative plasma [figure [2] (1)]. 
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Figure 4. PIC simulation results (He, 1 atm, 13.56 MHz, 2 mm gap, 330 V): Spatio- 
temporal plots of the (a) electron density, (b) space charge density, (c) electric field, 
(d) electron conduction current density, (e) electron heating, and (f) ionization rate, 
(g) shows the electric field profile at the time of max. ionization [vertical dashed lines 
in (a) - (f)] and the second term (drift field) of eq. ([3]). (h) shows the conduction 
current density in the discharge center and the applied voltage as a function of time. 
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Decreasing the driving voltage amplitude from 500 V to 330 V in case of atmospheric 
pressure micro discharges operated in helium under otherwise identical discharge 
conditions as shown in figure Q] significantly affects the electron heating and ionization 
dynamics, as shown in figure HI The electron density is approximately one order of 
magnitude lower (a) and the electric field (c), the electron conduction current density 
(d), the electron heating rate (e), and the ionization rate (f) are maximum in the bulk 
at a later phase within the RF period. This is caused by a higher bulk resistance due to 
the lower plasma density, which now dominates compared to the sheath impedance and 
reduces the phase shift between current and voltage to 3.6°, i.e. the discharge becomes 
more resistive pTTlH2]. Since current and voltage are almost in phase, maximum bulk 
electric field, current, heating, and ionization are observed at a later phase within the 
RF period. The electron heating is again completely dominated by ohmic bulk heating. 
The heating and ionization rates peak in the discharge center and no longer at the 
sheath edges, since the conduction current is only high in the centre. The high electric 
field is still caused by a low DC conductivity due to the high collision frequency, i.e. 
the second term of equation ([3]) reproduces the field obtained from the simulation at 
the time of maximum ionization well. Again, the slope of the electric field profile in the 
bulk is correlated with a small local space charge [40J. 

4. Conclusions 

Electron heating and ionization dynamics in capacitively coupled atmospheric pressure 
microplasmas operated in helium at 13.56 MHz and different voltage amplitudes were 
investigated by PIC simulations and semi- analytical modeling. The results were 
compared to electropositive argon and electronegative CF 4 macroscopic low pressure 
capacitive RF discharges. The electron heating dynamics in atmospheric pressure 
microplasmas is found to be dominated by the ohmic bulk heating of electrons due 
to high electric fields in the discharge center at the phases of maximum current in the 
RF period. The model reveals the physical origin of this f2-mode by identifying the 
high electric fields in the discharge center to be predominantly drift fields caused by 
a low conductivity due to the high electron-neutral collision frequency at atmospheric 
pressure. This heating mode is similar to the Drift-Ambipolar heating of electrons in 
low pressure electronegative macroscopic capacitive discharges operated in CF4, where 
a high bulk electric field is caused by a low conductivity due to a low electron density 
caused by the high electronegativity. 

In atmospheric pressure microplasmas, the heating and ionization dynamics are found to 
be affected by the amplitude of the driving voltage waveform. At high amplitudes, the 
plasma density is high and the bulk resistance is comparable to the sheath impedance, 
so that current and voltage are approximately 45° out of phase. Maximum ionization 
is observed adjacent to the sheath edges due to local maxima of the drift field caused 
by the ion density profile, that decreases towards the electrodes. These maxima of the 
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ionization are not caused by direct interaction of electrons with the oscillating sheath 
electric fields such as observed in low pressure macroscopic electropositive discharges. 
There is no cooling of electrons during sheath collapse. At lower driving voltage 
amplitudes, the plasma density decreases, so that the bulk resistance increases and 
the discharge becomes more resistive. Consequently, the phase shift between voltage 
and current decreases and maximum ionization is observed in the discharge center at a 
later phase within the RF period. 

These results might improve the understanding of measurements of the spatio-temporal 
emission in /z-APPJs [1~4T[T7] by Phase Resolved Optical Emission Spectroscopy p.8j and 
serve as a basis for a better understanding of chemical processes in such plasmas. Future 
investigations of the Jl-mode in electronegative atmospheric pressure microplasmas as 
well as experimental verifications of the observed dependence of the phase shift between 
current and voltage on the driving voltage amplitude are clearly required. 
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